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Abstract 

General substitution rules for non-periodic rhomb tilings are derived. From the 
requirement that all substitution tiles consist of a discrete number of prototiles, 
it follows that a substitution tile with angle S7r/n must be built out of pairs of 
prototiles with angles (s±t)7r/n, except if t = 0. In addition, we require that a 
discrete number of prototile edges must fit between the beginning and endpoint of 
the substitution tile edge. By comparing the total area of the discrete number of 
prototiles constituting the substitution tile and the total area derived from the as¬ 
sumed edge shape, a set of substitution rules is derived. The generalization involves 
the introduction of prototiles with a negative area or subtraction tiles. 


1 Introduction 

The interest in non-periodic tiling is shared by seemingly different groups of people, such 
as crystallographers, mathematicians, artists, or just individuals interested in the beauty 
of geometric patterns. The discovery of the first quasicrystals [1] in the eighties of the 
last century triggered many solid state physicists and chemists to explore solids having 
structures intermediate between crystallographically ordered and completely amorphous 
materials. Mathematicians designed new concepts to distinguish and classify structures 
showing aperiodic order. An extensive description of the achievements up to now is given 
in the excellent book by Michael Baake and Uwe Grimm |2]. 

A powerful way to generate an aperiodic tiling is by the repetitive substitution of a set 
of basic tiles, the socalled prototiles. The most commonly used tiling substitution rules 
contain a prescription of the way in which the linearly expanded prototiles are replaced 
by a number of original prototiles. This type of substitution rule is called geometric in[3], 
to distinguish it from a combinatorial one. The most famous substitution tilings based on 
pairs of rhomb tiles with opening angles svr/n, s = 1,2 are the Ammann-beenker (n = 4) 
[1] and the Penrose (n = 5) [3] Pj tiling. For the Socolar (u = 6) [7] and Goodman- 
Strauss(u = 7) [S] tiling a set of three rhomb tiles are needed. The Goodman-Strauss 
tiling was later generalised to arbitrary u > 3 by Harriss |9]. His general substitution 
tile is composed of 16 prototiles, having edges with either a dent or a dimple in the 
middle with an edge angle of 7i/n. Special substitution rules had to be designed for the 
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smaller substitution tiles. In the following section we will show that an even more general 
approach is possible. 


2 A general model to construct rhomb substitution 
tiles 

The model is based on the observation that the combined area of the pair of rhomb 
prototiles and is proportional to the area of the rhomb prototile T”, the propor¬ 
tionality factor 2cos(t7r/n) being independent of s. Consequently, substitution tiles 
may be constructed from a combination of prototiles T” and pairs of tiles If no and 
Ut are their numbers respectively, the area of the substitution tile is 

S = no + Tjt2nt cos(t7r/n) (1) 



Figure 1: Substitution tile edge structure. The edge angles Ui occur in ± pairs or are zero. 
In this example we chose the edge angles a* to be 0, —27r/9 , 7r/9, 27r/9, —7r/9 respectively. All 
edges have congruent shapes. The lower and upper left edges are related by a rotation over the 
opening angle s7r/9 with respect to the left corner. Similarly, the lower and upper right edges 
are related by a rotation over the opening angle with respect to the right corner. Opposite edges 
are related by a translation. 

A second requirement for a tiling of the entire plane is to realize proper edge sub¬ 
stitutions. We will assume that all four substitution tile edges have the same shape. 
Neighboring edges at the opening angle are related by a rotation over that angle and 
opposite edges are related by a translation (£gure[^. This edge arrangement also ensures 
that the substitution tile area is equal to the inflated rhomb area, and, therefore, S is 
equal to the areal scaling factor. The angles between the outer prototile edges and the 
substitution tile rhomb edge will be called the edge angles ai. For now, we will assume 
that overhangs are not allowed and \ai\ < 7r/2. The length of the substitution tile rhomb 
edge 


L = Tji cos ai 


( 2 ) 


is the inflation factor of the rhomb tiles. Because the areal scaling factor S is the 
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( 3 ) 


square of L, equations and can be combined into 

77-0 + cos{t7i/n) = -SjSj{cos(Q!j + aj) + cos(aj — aj)} 

This equality can only be satisfied, if the arguments a* + aj and — aj are both equal 
to an integer times 7i/n for all i and j. There are two solutions: either all angles a* are 
equal to an integer times tt/u, or all of them are equal to a half-integer times tt/u. 

Because the beginning and end of the substitution edge have to be at the endpoints 
of the substitution tile rhomb edge, the following relationship between the edge angles a* 
should be met: 

EjSmQfj = 0 (4) 

A general solution is that the edge angles occur in ±-pairs or are zero , and 
equation becomes 


no -l- 2T,tnt cos(t7r/n) = {mo -l- 2'Ermr cos(r7r/n)}^ 


(5) 


or ^ 

no-|-2Stnt cos(t7r/n) = {2Srm^_,_i cos((r-|- -)vr/n)}^ (6) 

Equations or determine the type and number of prototiles from which the 
substitution tile can be constructed, once the shape of the substitution tile edge has been 
chosen. In view of the above considerations, this edge shape may be characterized by 
a sequence of integers or halhntegers, the edge sequence (ki)^, dehned by ccj = kiTi/n, 
-n/2 <ki<n/2^5\. 

If the hnite edge angles are present as pairs in accordance with equation always a 
valid solution for the substitution tile is obtained, because both sides may be written as 
a sum of cosine terms having even valued coefficients. The pairing of the edge angles, 
therefore, guarantees that the substitution tiles are composed of an integer number of 
prototiles. 



Figure 2: General substitution scheme for a substitution tile S'” with a given edge shape. The 
prototile at position i,j is Tg+ki-kj, where kiir/n and kjTr/n are the edge angles at the upper 
and lower left edges respectively. 

^■There are also special solutions. For instance, if one requires that the sum of three terms are zero, 
one finds that ai + a 2 = tt/S and ai + = —tt/S. This solution is valid if n is a multiple of 3. An 

example satisfying this condition is the Lord tiling, having edge angles Oi = 02 = tt/G and 03 = — 7r/2 
[H]. n = 3 in this case, and the edge sequence is (|, — |, |). In this paper, however, we will only consider 
the more general ± pairing condition. 
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Equations]^ or [^constitute a connection between the prototile edge angle pairs rrik and 
the numbers of prototiles in a substitution tile rit, not their arrangement. The relations 
do not guarantee that a consistent set of substitution tiles can be found. However, in the 
following we will show that a general set of substitution rhomb tiles can be constructed 
for arbitrary n and for an arbitrary substitution tile edge shape. 

We start with a construction of the circumpherence of the tile S'” as described earlier 
and illustrated in figure [Tj Next, copies of the edges are translated to the breaks of 
neighboring edges. If the breaks of the upper and lower left edge are indexed as i G (1, 2..n) 
and j G (l,2..n) respectively, starting at the left corner as indicated in figure j^, one 
obtaines a grid of vertices (i, j), at which four prototiles meet. The one bounded by the 
vertices (^, j), {i + 1, j), {i,j + 1) and {i + 1, j + 1) is a prototile of the type The 

vertices at diagonal positions are occupied by tiles T” , whereas one can End pairs of tiles 
T'^±[ki-k-) off-diagonal positions (i,j) and (j, i). This general substitution rule may be 
represented by the following matrix 



n~'n 
^ s 

n~'n 

s+ki—ko 

'~pn 

s+k2-kQ 

n~'n 

S + kn-l- 

-ko 


n~'n 

s-\-ko-ki 

n~^n 

s 

n~'n 

-^s+k2-ki 

n~'n 

S + kn-i- 

-ki 

5r = 

n~in 

s-\-ko-k2 

'Y'n 

s+ki—k2 

'~rn 
^ s 

n~'n 

S^kn-l- 

-^2 


'Y'n 

_ s-\-ko-kn-i 

'Y'n 

S+fci-fcn-l 

n~'n 

S-\-k2-kn-l 

n~'n 

• * s 



For later use one should note, that the prototiles parallel to the substitution edges, 
i.e. the rows or columns of the matrix, form worms, and the edges of the worms have 
shapes identical to the edge shape of the substitution tile. 

The prototiles are allowed to have indices s ± t < 0 or > n. These prototiles will have 
negative areas, meaning that they have to be subtracted from the tiling. We consider a 
tiling of the plane to be a legitimate one, if in the end there are no holes or overlaps. 
So, negative or subtraction tiles are allowed, if they remove all overlaps between tiles and 
do not leave holes in the tiling. In one of the next sections, we will reason, that this is 
presumably the case for substitution edges without loops. Also the zero area prototiles for 
which s ± t = 0 or n play a important role in our scheme and cannot simply be neglected. 


3 Substitution tiles with mi = 1 


We will now discuss one of the simplest examples of this general scheme, i.e. the one with 
mi = 1 and m>i = 0. From it follows that 

no = ml + 2ml, = 2momi, n2 = ml (8) 


For mo = 0 we get no = 2 and ni = 1 and the 2x2 general substitution tile is 


S 


n 

s 




■s-2 


n~>rL 

-^s+2 



(9) 


The edge sequence is (1, -1) and the inflation factor L = 2cos(7r/n). 

In figure all first generation substitution tiles, constructed using the general scheme 
outlined in the previous section, are shown for n = 9. Negative areas are colored white 
in a gray background. The s = 0, 1, 8 and 9 substitution tiles are constructed in exactly 
the same way as the general tile, i.e. the s tiles are rotated over ^nir/n and connected 
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Figure 3: mi = 1 tiles for n = 9. The s = 0,1,8 and 9 tiles have been constructed using 
negative area prototiles. 

with the unrotated s ± 2 tiles. For s = 0 two zero area prototiles Tq, which are lines of 
length 2, a prototile and a T ”2 subtraction prototile have to be used. The normal and 
subtraction tile annihilate each other, and the substitution tile is reduced to a zigzag line. 
This is most easily seen by considering tiles T^, and T^j -2 letting S approach 0. For s = n, 
on the other hand, the T ”^2 the T ^_2 tiles do not overlap and the substitution tile 
consists of a positive and negative part. However, the total area is zero, as it should. For 
s = 1 the negative area T"^ tile cuts pieces off the other three prototiles. For s = n — 1 
the two prototiles T^_i partly overlap. This overlap is compensated by the negative area 
tile, leaving a fractional negative region shown in white. An essential result is that 
all these tiles have the same dents and dimples as the general tile and can be used to 
construct larger substitution tiles. 



Figure 4: mi = 1 tiles for n = 5. White parts have a negative value and are to be subtracted 
from the tiling. Four different sets of substitution tiles can be constructed for the second 
generation. 

The prototiles T” and have identical rhomb shapes and, therefore, one obtains 
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« = 8, s odd 
s = -\ 1 

r‘ generation substitution tiles 







2'“‘ generation substitution tiles 




Figure 5: mi = 1 tiles for n = 8. Odd and even s tiles form separate tile subsets, which do not 

mix.The even s subset is identical to the mi, n = 4 tile set. 
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an identical object if one replaces the other. The same holds if a prototile is rotated 
over 71. However, this is different for higher generation substitution tiles. Although the 
substitution tiles and S'”_g have the same basic rhomb shapes, their edge conhgurations 
are different. A replacement or a twofold rotation of one or more of the four prototiles, 
leads to 16 possible tiles with different edge shapes. The reasoning is as follows: if we £11 
up a next generation S'" tile, there are four ways to put a T" or T^_s tile at the (1, 1) 
position. Subsequently, two out of four fit at the (1, 2) and (2, 1) 

position. Finally, only one out of four T" or type tiles matches the edges already 
present at the (2, 2) position. The total number of tile shapes is, therefore, sixteen. 
However, only four of them are consistent with the requirement that the four edge shapes 
are the same and related to each other by a translation or rotation as defined in figure [T] 
These are: a) the basic substitution tile, b) a substitution tile in which all constituent tiles 
in a) have been rotated over vr, c) a substitution tile in which the s ± 2 tiles are replaced 
by n — s =F 2 tiles and d) a substitution tile in which all constituent tiles in c) have been 
rotated over tt. Another, simpler way to arrive at the four substitution rules is to consider 
the edge sequence of a tile. There are two ways to connect an edge between two corners of 
a rhomb. Consequently, the number of different edge shapes for a substitution tile is 2^, 
if N is the number of prototile edges along the substitution tile edge, which is two in the 
case considered here. The four second generation substitution tile varieties are shown for 
n = 5 in figure and n = 8 in figure The first two are mirror symmetric, whereas the 
last two are asymmetric, but mirror images of each other. All four sets contain tiles with 
negative area parts. These negative areas account for overlapping dimples of neighbouring 
edges. They will naturally disappear if the next larger substitution tile is constructed, 
i.e. overlapping dimples will be filled up with dents of neighboring tiles. In some cases, 
notably the s = n and s = n — 1 tiles, new negative areas will appear at one end of the 
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Figure 6: Edge shapes for three generations of n = 5 2 x 2 full dent substitution tiles. 



Figure 8: Aperiodic n = 5 von Koch Rhomb Tiles. 


tile. Also in most of the basic substitution tiles of set a), the upper part is a negative area 
part. Second generation substitution tiles with 0<s<n — Ido not have negative area 
parts for sets b), c) and d). In set b) this is due to the fact, that the s = n and s = n — 1 
constituent tiles are used in a rotated position, so that no negative area parts stick ont. 
In sets c) and d) the s = n and s = n — 1 prototiles have been replaced by the s = 0 and 
s = 1 prototiles respectively, not having negative area parts by themselves. 

Because the s-values of the prototiles in substitution tiles belonging to set a) and b) 
differ by 2, the even and odd s snbstitntion tiles constitnte separate snbsets, both fully 
tiling the plane. For even n this is also trne for sets c) and d), whereas for odd n the even 
and odd tiles intermix. 
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Some of the prototiles are cut into smaller and rather complicated shapes by the 
subtraction procedure. The latter is not unique, because the hnal surface is obtained by 
cutting the shape of the negative area out of one of two overlapping dents. In hgure 
this region has been hlled in quite arbitrarily by showing a linepattern consistent with 
the symmetry of the tile. Although we cannot prove it rigorously, it appears that it will 
be possible to replace the mutilated tiles by half-prototiles (triangles) for n=5. 

Each new generation of mi = 1 substitution tiles will be fourfold in the same way 
as was described for the second generation. For each new generation, an edge segment 
is replaced by a two-segment dent or dimpel. Already after a few generations a large 
number of different edge shapes are generated. In hgure the edge shape is shown for 
three generations for n = 5. Sometimes, overhangs or loops show up. The latter are due 
to overlapping neighboring dimples causing negative area parts at the circumpherence of 
a substitution tile. In hgure this is seen to occur once. The way a next generation 
edge shape is obtained is reminiscent of the construction of a von Koch curve uni. An 
important diherence is that in the original von Koch construction each segment is replaced 
by four instead of two new segments. The actual von Koch curve is obtained if substitution 
rule b) of hgure is applied for each new generation. Eight generations of these von Koch 
edges are shown in hgure and the corresponding von Koch tiles in hgure A complete 
tiling of the plane may be generated using either the odd or even s tiles. 


4 Substitution tiles with mi = 1 

2 

In section!^ we found that there are two types of edge substitution rules, one with integer 
and one with half-integer edge angle fractions of tt, eqs. and respectively. In this 
section we will discuss the second case. From equation]^ 

no = 2m|; rii = m\] n. i = 0 (10) 

2 2 2 



Figure 9: mi = 1 substitution tiles for n = 9. 

2 

If mi = 1, it follows that no = 2 and ni = 1. Thus, like in the case of mi = 1, four 
2 

prototiles are involved, but the non-congruent tiles have s-values dihering only by one 
instead of two, or 


S 


n 

s 


rj-in 


• s —1 


s+1 


The edge sequence is (|, — |) and the inhation factor L 


cos 7r/2n. 


( 11 ) 









s=-4 -3 -2 -1 0 1 2 3 4 5 

0 0 o o — o ♦ ♦ 0 I 



2”‘‘ generation m„ 2 =l substitution tiles 


Figure 10: mi=ln = 5 (Penrose Rhomb) tiles. White denotes a negative area. Like in the 
2 

772-1 = 1 case, four different 2nd generation tile sets may be constructed, two mirror symmetric 
sets a) and b) and two asymmetric sets c) and d). The edge sequence of the second symmetric 
set (b) is identical to the one of the Harriss re = 5 set of tiles. The odd s subset of c) and d) is 
the so-called binary tiling [8] 


In figure]^ the 1st generation substitution tiles for n = 9 are shown. In contrast to the 
mi = 1 substitution tiles (£gure[^, negative area tiles are not needed for the construction 
of positive area tiles. Also in this case, there are four possible edge shapes for the second 
generation of substitution tiles, depending on how the first generation substitution tiles 
are combined (see hgure 10 for mi = 1, re = 5 tilings). The s = re — 1 substitution tile 
of set (a) contains a negative area, because the dimples at neighboring edges overlap. In 
higher generation substitution tiles this negative area will also appear in the tiles with 
smaller s-values. For set b) the negative half of the s = n constituent tile does not show 
up, because it removes the overlap of two s = re — 1 tiles. Also in the higher generation 
tiles negative area will not be visible. In fact, the shape of the edge of the tiles of set b) 
is identical to the shape of the Harriss tiles. The tiles of set c) and d) are, like in the 
mi = 1 case, mirror images of each other. For odd re, they can be split up into two sets, 
both giving a full aperiodic tiling of the plane, one with odd and one with even s-values. 
The one with even s tiles is special because no s = re constituent tiles are needed and, 
therefore, no negative tiles are involved. By neglecting the the s = 0 tiles, the set of 
prototiles is reduced to a conventional one. For re = 5 the tiling based on the even s set, 
i.e. s = 2 and s = 4 is already known as the ’’binary” or Langon Billard (LB) tiling |llj . 
The odd s set contains the s = n tile. Its role is to remove a patch of prototiles at one 
side of an s = 1 tile and add it to the other side. A tiling results which is very similar to 
the LB tiling. The prototile set for even re comprises all s tiles from 0 to re. Although the 
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negative part of the s = n tile is involved in the substitntion rnle, the prototiles are not 
cut up and a tiling consisting of complete prototiles results. 

The mi = 1 substitution tiles are identical to a subset of rrii = 1 tiles with a doubled 


n-value and with even s-values. So, the tiles shown in figure 10 are also obtained as an 


even s subset of the mi = 1, n = 10 case and the even s subset of the mi = 1, n = 8 tiles 
shown in figureis identical to the mi = 1, n = 4 tile set. 

5 Substitution Matrices. 

In this section we want to reformulate the rhomb substitution model in terms of the edge 
and tile substitution matrices. If overhangs are included, the prototile edges in a tiling 
will point into 2n directions. Each of these is replaced by a number of prototile edges in 
orientations determined by the edge sequence, i.e. uiq in the same, m„ in the opposite 
direction and m^ with k G (1,2, ..,n — 1) in directions differing by +/c7r/n and —kir/n. 
The edge substitution matrix, therefore, is 


M = 


■ mo 

mi 

m2 

■ ■ ■ mn-l 

mn 

... ms 

m2 

mi 

mi 

mo 

mi 

■ ■ ■ mn-2 

mn-l 

m4 

ms 

m2 

m2 

mi 

mo 

■ ■ ■ mn-3 

mn -2 

... ms 

m4 

ms 

rrin-i 

mn-2 

mn-3 

... mo 

mi 

. . . mn -2 

mn-l 

mn 

m„ 

mn-i 

mn-2 

mi 

mo 

■ ■ ■ mn-3 

mn -2 

mn-l 

ms 

1714 

ms 

• • • mn -2 

mn-s 

... mo 

mi 

m2 

m2 

ms 

1714 

■ ■ ■ mn-l 

mn -2 

mi 

mo 

mi 

mi 

m2 

ms 

... ITln 

mn-l 

m2 

mi 

mo . 


( 12 ) 


A tile with index s is substituted by no prototiles with index s, n„ with index s + n 
and Ut with index s + t and s — t, with t G (1,..., n — 1) and s G (0,..., 2n — 1). So the 
substitution matrix S is 


S = 


■ no 

Til 

n2 

■ ■ ■ nn—i 

nn 

... ns 

n2 

ni 

ni 

no 

ni 

■ ■ ■ nn-2 

nn—1 

n 4 

ns 

n2 

n2 

ni 

no 

■ ■ ■ nn —3 

nn-2 

... ns 

TI4 

ns 

'^n—l 

nn—2 

nn —3 

... no 

ni 

• • • nn—2 

nn—1 

nn 

Un 

nn—1 

nn—2 

... ni 

no 

■ ■ ■ nn—s 

nn—2 

nn—1 

ns 

714 

ns 

• ■ ■ nn—2 

CO 

1 

... no 

ni 

n2 

n2 

ns 

714 

■ ■ ■ nn—1 

nn—2 

... Til 

no 

ni 

. ni 

n2 

ns 

... Tin 

nn-1 

... ns 

ni 

no . 


The relation between the tile and edge substitution matrices is 

S = 


(13) 


(14) 
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This matrix equation may be used to calculate the numbers of prototiles in a substi¬ 
tution tile for a given edge shape instead of equations . The Uj are given by the product 
of the hrst row and the j-th column 




— Sij+i — j G (0, 


( 16 ) 


Both M and S are circulant 2n x 2n matrices na. Consequently, a shorthand notation 

is 


of equations 12 and 13 


M = circ{mo,mi,m2, m 3 , m 2 , mi) 


(16) 


S = circ{no,ni,n 2 , ...,n„_i,n„,ns, ^ 2 , ui) (17) 

All circulant matrices are known to have the same set of normalized eigenvectors 

1 


vi = 


y(2n) 


-( 1 . 


,(2n- 




with e = expijii/n) and I G (0,1,..., 2n — 1). The eigenvalues of M are 

\j = = fno + 2S^“Jmfc cos(7rjfc/n) -|- (—l)-^m^ 


(18) 


(19) 


, and because of relation 14, those of S are A^. The eigenvector Ai is equal to the inflation 
factor 


L = mo - rUn + (m^ - mn-k)cos{k'K/n) (20) 

A substitution tiling can only be a model set for a quasi crystal if its inflation factor is a 
Pisot- or PV-number [13], because a model set is point diffractive na. L is a PV-number, 
if the absolute value of all its conjugates is less than 1. The conjugate eigenvalues are the 
ones for j coprime to n. Using the above formulae we hnd that the inflation factors are 
PV-numbers in the following mi = 1 or m 2 = 1 or socalled single dent cases: 


mo 

mi 

m 2 

n 

0 

1 

0 

5 

1 

1 

0 

4,5,6,9 

2 

1 

0 

4,6 (Harriss) 

1 

0 

1 

5,7 

2 

0 

1 

5 


Table 1: Single dent substitution tiles with PV inflation factor. 


The edge substitution matrix for a halhnteger edge sequence can be obtained by dou¬ 
bling the n value. The fractional indices have to be doubled as well and become the m^ 
values for odd k, whereas the m^ for even k are zero. From table [T] it is clear that the half 
integer single dent substitution tiles will not have PV inflation factors. 
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6 Relationship with known rhomb tilings 

In all known rhomb tilings only positive area tiles are used. Nevertheless, some of the rules 
derived in this paper can also be applied. The examples we will discuss are the Ammann 
Beenker(AB) tiling, the Penrose (P) tiling, the binary or Lancon Billard (LB) tiling and 
the Harriss n-fold generalization of the Goodmann-Strauss tiling. The inflation factors of 
many known rhomb substitution tilings can be calculated using either formula or 
One requirement is that the non-zero prototile edge angles occur in ± pairs, another that 
the area of the sum of the prototiles in the substition tile is equal to the substitution 
rhomb area. Than, the number of prototiles per substitution tile and, therefore, the 
substitution matrix can be calculated using equation by neglecting the zero area tiles, 
and subtracting the negative tiles from the corresponding positive ones, or = ns—n 2 n-s- 
This set of positive prototiles can be made even smaller by adding the ones having identical 
basic rhomb shapes or n" = Ug + Un-s — nn+s — n 2 n-s- 

The most famous examples are the AB (n = 4) and P(?7, = 5) tilings in which the full 
set of 8, respectively 10 tiles have been reduced to a set of only 2 prototiles. The AB 
substitution tile has a (0, 1, -1) edge sequence. However, the edge conhguration is different 
from the one in £g. and differs for both prototiles. The P substitution tiles have two 
different pairs of edges, a couple with a (1, -1) and a couple with a (0, 2, -2) edge sequence. 
These edges have the same length, because 2cos(7r/5) is equal to 1 -|- 2cos(27r/5). The 
corresponding sets of m-values, (mi = 1), and (mo = l,m2 = 1) respectively, also yield 
the same substitution matrix S for n = 5. 

Some of the known rhomb tilings are in a way special cases of our general scheme. 
For instance, the binary or Langon Billard(LB) tiling. The LB substitution tiles have 
an mi = 1 edge. As was discussed in section the even s tiles of the n = 5 case are 
identical to the substitution tiles of the LB tiling, if substitution rule c) or d) is followed. 
Our treatment shows, that the LB tiling can be generalized to odd n values different from 
hve. The prototile set consists of all even s tiles. The set of odd s tiles produces another 
very similar tiling. Although the s = n tile with a positive and negative part is involved, 
the prototiles in the substitution tiles are not cut up. For even n, the full set of prototiles, 
odd and even, are to be used to construct an LB-like tiling. 

Another special case of our tiling scheme is the Harriss tiling for arbitrary n > 3 of 
which the Goodmann-Strauss tiling (n = 7) is an example. The edge sequence of the 
Harriss tiles is (0, 1, -1, 0), i.e. mo = 2 and mi = 1. The corresponding rit values are 6, 
4 and 1 for f = 0, 1 and 2 respectively. By replacing the s ± 1 tiles by n — s =F 1 tiles, as 
Harriss did, the set of prototiles is reduced to the even s prototiles in the odd n case. From 
the perspective of our scheme it has the advantage that the s = n tile with a negative 
region is replaced by the s = 0 tile whithout a negative region. For even n, Harriss had 
to introduce a number of additional tiles having extra dents or dimples. (Note, that the 
index p used by Harris is equal to n — s) 

Finally, the tiling with 11-fold rotation symmetry found by Maloney with computer 
assistence, is based on a set of only five complex prototiles with a (l,-l,-3,3,0,2,-2,-l,l,0,- 
5,5,-3,3,-1,l,4,-4,2,-2,0,-1,l,2,-2,-3,3,0,4,-4,-l,l,2,-2,0) edge sequence [15]. The correspond¬ 
ing set of (mt)Q is (5,5,4,3,2,1). The positive and even s tiles, are apparently sufficient to 
tile the plane. 
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7 Gaps and Overlaps 

A consequence of the admittance of negative prototiles in our general model is that mem¬ 
bers of the substitution tile set might contain negative parts, gaps and overlaps. In the 
following we will analyze the conditions under which these defects disappear in tilings. 





a) b) c) 


Figure 11: Tile construction in terms of rows of prototiles or worms, a) Row of prototiles 
along an edge. At the crossing point of the upper and lower worm edges the exterior negative 
and positive prototile parts compensate each other and the sign of the interior prototile (parts) 
changes, b) By adding the next row of prototiles one gets a double worm following the same 
rule, i.e. the interior parts are either negative or positive and the exterior is empty, c) Once the 
tile is completed, it will consists of positive and negative regions, if its edges cross. The example 
shown here is a T® tile with a (1,2,3,-3,-2,-1) edge sequence. 


As was suggested earlier, the general substitution rule for an iV x V substitution tile 
may be formulated in terms of worms of N prototiles. The procedure is illustrated in 
hgure n The worms are bounded by two parallel copies of a substitution tile edge and 
separated by one edge unit along neighboring edges. The interior of the worm, i.e. the 
space between the two edges, is either hlled with positive or negative tiles or tile pieces, 
provided the edge does not contain loops. The sign changes at the crossing points of 
the two. Adding a neighboring worm, one gets a double worm with congruent edges, 
and interior regions which are either positive or negative. If the tile edges do not cross 
anywhere after completion of the tile, i.e. if the circumpherence does not contain loops, 
the hlling of the tile is a single positive layer of tiles or tile pieces. If they do cross 
somewhere, there will be regions hlled with a layer of negative tiles or tile pieces. To fully 
tile the plane we may repeatedly inhate a positive tile and apply the edge substitution rule. 
This involves the substitution of the inhated edge segments by the initial edge sequence. 
Since the edge substitution may be done in two orientations related by a twofold rotation, 
each time 2^ diherent next generation edge sequences are possible. The new extended 
edge contains segments, g being the substitution generation. Consequently, the worms 
consisting of prototiles are the new building blocks of the tile. If we do not admit 

loops in the edge sequence, the interior of the multiply inhated tile will consist of positive 
and negative regions separated by the crossing points of the four edges, but no gaps or 
overlaps. 


8 Summary 

A general scheme for constructing a non-periodic rhomb tiling was derived for arbitrary 
n > 3. It admits zero and negative area tiles to ensure generality. If the edge sequence 
does not contain loops, such a tiling will consist of single layered positive and possibly 
negative regions, but there will be no gaps or overlaps in the tiling. Overlaps occurring in 
the tiling process are always annihilated by negative or subtraction tiles. This generally, 
but not always, leads to the cutting up of tiles. 
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The corresponding edge substitution matrix can be expressed as a circulant matrix 
M, the inflation factor being one of its eigenvalues. The tile substitution matrix S is 
equal to the square of M. For a given edge sequence this equality determines the type 
and number of prototiles in a substitution tile. A possible arrangement of the prototiles 
is obtained by drawing a grid using parallel copies of the edge. This arrangement can be 
expressed in the form of a matrix. 

The smallest (2 x 2) substitution tiles were discussed in detail. Two different types 
can be distingushed, those with an integer and those with a halflnteger edge sequence. 
The latter are of special interest, because some of the tilings do not contain tiles cut up 
by the application of subtraction tiles. An example is the LB tiling (n = 5), which can 
be generalized to arbitrary n. 
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